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Abstract 

We have searched for lepton flavor violating semileptonic t~ decays using a data sample of 153.8 
fb” 1 accumulated with the Belle detector at the KEKB e + e _ collider. For the six decay modes 
studied, the observed yield is compatible with the estimated background and the following upper 
limits are set at the 90% confidence level: B(t~ —> e~r ]) < 2.3 x 10~ 7 , B(t~ —> fx ~r j) < 1.5 x 10 -7 , 
B(t~ —> e _ 7r°) < 1.9 x 10~ 7 , B(t~ —> < 4.1 x 10~ 7 , B(t~ —> e~r]') < 10 x 10 -7 , and 

B(t~ —> n~r]') < 4.7 x 10 - '. These results are 10 to 64 times more restrictive than previous limits. 

PACS numbers: 11.30.-j, 12.60.-i, 13.35.Dx, 14.60.Fg 
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INTRODUCTION 


Processes with Lepton Flavor Violation (LFV) provide some of the most promising 
searches for physics beyond the Standard Model (SM). In the charged lepton sector LFV 
is usually considered for purely leptonic processes such as radiative decays of /i and r [l| 
or their decays into three charged leptons j2]. The semileptonic r decays r” — * t~ M° 
(where i~ = e”,yU” and M° = 7r°, 77 , 77 ') provide another good source of information about 
possible LFV mechanisms from supersymmetry (SUSY) motivated models 4| to models 
with heavy Dirac neutrinos jfjj, [ 6 ]; they also allow more general bounds to be placed on the 
scale of new physics |?]. Experiments at B factories, where r lepton pairs are copiously 
produced, have substantially increased the sensitivity to LFV decays, bringing it close to 
various theoretical expectations. Recently we performed a search for the decay r” —> n~rj 
based on a data sample of 84.3 fb ” 1 and placed a new upper limit on its branching fraction, 
£>(r” —> /U 77 ) < 3.4 x 10”' | 8 ], 28 times better than the previous limit from CLEO 9j and 
comparable to predictions for some SUSY models {3,0]. 

We present here an updated result for r” —> /i”^ and new searches for the decay modes 
r” —> /r”7r°, r” —> 77 ', r” —> e”7r°, r” —> e — 77 , and r” —* e~r)', based on a data sample of 

153.8 fb” 1 , equivalent to 137.2 xlO 6 t + t” pairs, collected at the T(4S') resonance with the 
Belle detector at the KEKB energy-asymmetric e + e collider na Due to this much higher 
integrated luminosity our sensitivity to the decay modes r” —> £~n 0 and £~r/ is considerably 


better than in previous searches with MARK II |ll| . Crystal Ball 

J 


13, ARGUS (13| and 

Unless otherwise 


CLEO J9f. The decay modes r” —> £~r]' are studied here for the first time 
stated, charge conjugate modes are implied throughout the paper. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic held. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons 
(KLM). A more detailed description of the detector can be found in Ref. Q. 

For Monte Carlo (MC) studies, the following programs have been used to generate back¬ 
ground events: KORALB/TAUOLA 15| for t + t” processes, QQ J16] for BB and continuum, 
BHLUMI Q for Bhabha, KKMC [3] for /i + /i” and AAFH jiff for two-photon processes. 
We use the following sizes of MC samples: 428 fb ” 1 of generic r + r”, 232 fb ” 1 of qq with 
q = u,d,s, and 160 fb ” 1 with q = c, 90 fb ” 1 of BB , 119 fb ” 1 of 9.1 fb ” 1 of Bhabha 

and 205 fb ” 1 of two-photon processes. Signal MC is generated by KORALB/TAUOLA. 
To take a range of possible r” —> £~M° production mechanisms into account, we generate 
events under three different assumptions: uniform angular distribution in the r rest frame, 
V — A interactions, and V + A interactions. We initially assume all r decays to have a 
uniform angular distribution and take the other hypotheses into account later. The Belle 
detector response is simulated by a GEANT3 |13 based program. 


DATA ANALYSIS 

We follow the same principles of event selection as those in the r” —> ii~rj search |sj]. 
Kinematical variables with a CM superscript are calculated in the center-of-mass frame, 
and other variables are calculated in the laboratory frame unless otherwise stated. We look 
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for an event with the following particles: 

{t~ —» t~ + M 0 } + {r + —> (track)" 1 " + 777 + X} , 

where the system X is unobserved. In other words, the events should be consistent with a 
t + t~ event, in which the t~ decays into a lepton £“ and a pseudoscalar meson M° (signal 
side) and the r + decays into a charged track, any number n > 0 of photons, and one or more 
neutrinos (tag side). The charged track on the tag side should not be an electron (muon) 
if the lepton on the signal side is an electron (muon), in order to avoid contamination by 
Bhabha (/i + /i~) events. We reconstruct mesons in the following modes: tt° —> 77, 7 —> 77, 
77 —> 7 r + 7 r - 7 r°, rf —> 7r + n~r] 1 7 —»■ 77. When pseudoscalar mesons M° are reconstructed from 
77, the event has a 1-1 prong configuration. When the 77 is reconstructed from h + ti~h 0 , or 
the 77 ' is reconstructed in 7 / —> 7 r" l ' 7 r — 77 , 77 —> 77, the event has a 1-3 prong configuration. 
The signal side thus contains at least two photons in all cases. 

For the 1-1 prong selection, the candidate should contain exactly two oppositely charged 
tracks and two or more photons, two of which form a 7 r° or an 77 . We require the momentum 
of each track, p, and the energy of each photon, E 1 , to satisfy p > 0.1 GeV/c and 7% > 0.1 
GeV. The tracks and photons are required to be detected within the CDC acceptance, 
—0.866 < cos 6 < 0.956. In order to exclude Bhabha, p + p~ and two-photon events, which 
otherwise contribute a large background, the total visible energy in the CM frame is required 
to satisfy 5 GeV < E™^ < 10 GeV. The event is subdivided into two hemispheres by a 
plane perpendicular to the thrust axis evaluated from all tracks and photons satisfying the 
above requirements. 

The lepton flavor is identified based on likelihood ratios calculated from the response 
of various subsystems of the detector. For electron identification, the likelihood ratio is 
defined as V e = C e /(C e + C x ), where C e and C x are the likelihoods for electron and other 
flavor hypotheses, respectively, determined using the matching between the position of the 
charged track trajectory and the cluster position in the ECL, the ratio of the shower energy 
in the ECL to the momentum measured by the CDC, the shower shape of the cluster in the 
ECL, specific ionisation ( dE/dx ) in the CDC and the light yield in the ACC (2ll ]. For muon 
identification, the likelihood ratio is 7% = + Ck + C n ), where C w , and Ck are 

the likelihoods for the muon, pion and kaon hypotheses respectively, based on the matching 
quality and penetration depth of associated hits in the KLM 22|. For electron candidates 
Ve > 0.9 is required, identifying electrons with an efficiency r] e = (92.4 ± 0.4)%; the pion 
fake rate (the probability that a pion is misidentihed as an electron) n e = (0.25 ± 0.02)%. 
For muons 77 ^ = (87.5 ± 0.3)% and = (1.4 ± 0.1)% are obtained. If the track satisfies 
V e > 0.6 or Vfj, > 0.6 on the tag side, i.e. is lepton-like, the constraints 7 ?- 7 ag < 2 and 
0.7 GeV/c are applied to suppress events including photons from initial state radiation or 
beam background. Otherwise, no constraints are applied on the tag side. The lepton flavor 
identification on the signal side is postponed until the last stage of selection. 

The momentum of a 7 r° or 77 meson produced in a two-body r-decay is on average higher 
than that of n 0 /r] mesons from other sources, so that photons from 7T° or 77 are required to 
have a rather high energy E 1 > 0.30 (0.22) GeV in the case of £“ = To further 

reduce the backgrounds, the cosine of the opening angle between and 77 on the signal side 
must satisfy 0.5 < cos 6%_ 77 < 0.95. When reconstructing an 77 meson candidate, a 7 r° veto is 
applied: defining the resolution-normalized 7 r° mass S— 0.135 GeV/c 2 )/cr 7 y, we 
reject events where a signal-side photon 7 and a tag-side photon 7' satisfy —5 < S*°, < +5; 
the resolution a 77 , is in the range 5-8 MeV/c 2 . This 7 r° veto rejects 86 % of the BG events 
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FIG. 1: cos#™,,. miss distribution for the r —> p q decay mode. The data (points with error 
bars), total background MC (open histogram), non-r + r~ background MC (hatched histogram), 
and signal MC (shaded histogram) are shown. The selected region is indicated by the arrow. 


while retaining 75% of the signal. To ensure that the missing particles are neutrinos rather 
than 7 ’s or charged particles that fall outside the detector acceptance, we require that 
the direction of the missing momentum satisfy the condition —0.866 < cos < 0.956. 
Since neutrinos are emitted only on the tag side, the direction of the missing momentum 
Pmiss should be contained on the tag side. We define the angle ^thrust-miss between the 
thrust axis of the event (in the signal hemisphere) and the missing momentum vector, and 
require cos ^[hrust_ miss < —0.55 as shown in Fig. [I] At this level of selection the dominant 
background is that of t + t~ events (95%) whereas qq events and other sources constitute 
only 5%. The correlation between the missing momentum, p miss , and the missing mass 
squared, iss , is used to remove the generic r + r~ and qq continuum contributions: p m i SS > 
—2.616m% ss — 0.191 and p m i SS > 1.0m% ss — 1, where p m i SS is in GeV/c and rn miss is in GeV/c 2 . 

The yield of signal events is finally determined in the Mi-mo-AE plane, where Mi-mo 
is the invariant mass of the £~M° system and A E is the difference between the energy of 
the £~M° system and the beam energy in the CM frame. When deciding on our selection 
criteria, we blinded the signal region M r — 5ct 1 £^_ m0 < Mi - m o < M r + 5and —0.5 

GeV < A E <0.5 GeV. We define sideband regions of Mg_- M o and A E as M r — 10 <Jm ^_ m0 < 
Mi- M o < M r - 5^ m0 , M r + 5 < Mi- M 0 < Mr + IOct^^ and -10cr^ < A E < 


TABLE I: The resolutions of Mg - M o and A E for each mode. The “low” or “high” superscript 
indicates the lower or higher energy side of the peak; a m ( _ v;0 are in MeV/c 2 and o/\e are in MeV. 


Mode 

Subdecay mode 

low 

a M t -M0 

high 

^M ( - M o 

low 

a AE 

high 

a AE 

t~ — ► e~p 

V -> 77 

25.8 

15.3 

79.3 

34.7 


< 

< 

T 

o 

o 

1 

+ 

T 

c- 

17.5 

9.7 

44.2 

25.5 

t~ -*• p~q 

V “»• 77 

25.7 

14.2 

57.1 

35.2 


< 

< 

t 

O 

o 

1 

+ 

T 

Sr 

13.4 

9.6 

43.2 

22.1 

T~ —> e _ 7 T° 

C~ 

< 

T 

o 

25.7 

14.7 

69.5 

35.2 

■s 

1 

I 

1 

O 

< 

< 

T 

o 

23.6 

14.5 

69.8 

38.9 

r — > e rf 

?/ —4 7 T + 7 T p (q — > 77) 

22.0 

14.0 

62.7 

26.7 

t~ -*• p~rf 

rf — > TT + TT~p (7 — > 77) 

18.3 

11.7 

55.6 

25.9 
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TABLE II: Numbers of remaining events for each mode at the same selection stage as in Fig. 0 
prior to final cuts on the signal side (see the text). 


Mode 

Subdecay mode 

ArDATA 

iV side 


t~ —> e~rj 

V -* 77 

3 

1.47 ± 0.73 


7/ — ► 7T + 7r _ 7T 0 

5 

0.34 ± 0.11 

T~ -*• PW 

V * 77 

17 

14.4 ± 3.2 


r] -7- 7 r + 7 r _ 7 r 0 

7 

5.2 ± 3.6 

t~ —> e~n 0 

7T° —> 77 

2 

2.2 ± 0.9 

t~ —> fi~7 r° 

7T° —> 77 

22 

25.8 ± 4.5 

t~ —> e~ rj 

rf —> 7T + 7T _ 7 

12 

11.6 ± 2.3 

t~ -> pi~rf 

rf — > 7 r + 7 r ”7 

33 

22.9 ± 3.5 


—0.5 GeV, 0.5 GeV < A E < lOo^f? 1 , respectively. The resolutions, and cr^^ hlgh , 

of Mi ~ m o and A E, are evaluated from the distributions of signal MC around the peak using 
an asymmetric Gaussian shape to account for initial state radiation and ECL energy leakage 
for photons. In these expressions, M r is the central value of the reconstructed r mass for 
signal, evaluated in MC: it is consistent with the known r mass within 1.2 - 13.0 MeV for 
each mode. The resulting resolutions are summarized in Table |I| 

The resulting numbers of data and MC events in the M^-mo-AE sideband region after 
applying a loose mass requirement for n 0 /rj, —5 < S~ r / V < 5, are summarized in Table ITT! 
The Sfi and SJj distributions for the modes with a final state muon are shown in Fig 121(a) 
and (b), respectively. Within the statistical uncertainty the sideband data and background 
MC yields are consistent with each other in all four modes. 


At the last stage of the selection, the n° or rj candidate is required to satisfy —5 < SfJ' n < 
3, where S^° was defined above, and S ^ = (m 77 — 0.547 GeV/c 2 )/cx ^ 7 with cr^ 7 = 11-13 
MeV/c 2 . The conditions V e > 0.9 and p > 0.7 GeV/c are imposed on lepton tracks on 
the signal side. In addition, V e < 0.6 is required on the tag side for rejection of Bliabha 
or p + /i _ events. 

For rj —> 7 T + 7 r - 7 r° and rf —> 7i + 7r~ r r] decays with p —>• 77 , i.e. 1-3 prong modes, we search 
for events containing four charged tracks (net charge = 0) and two or more photons. Because 
of the higher multiplicity compared to the 1-1 prong modes, the detection efficiency becomes 
smaller. On the other hand, the additional reconstruction constraint in the 7/77' decay chain 
improves the background rejection power. The selection criteria are similar to those in the 
rj —> 77 case with the differences listed below. 

The minimum photon energy E 7 requirement is reduced from 0.1 GeV to 0.05 GeV, and 
the tighter conditions for photons from 7r °/r] are also removed, since the photons from these 
decay modes have a softer energy distribution compared to that in the 1-1 prong case. At 
least three tracks and two or more photons are required in the signal hemisphere. One track 
must be identified as an electron or a muon (' P e / M > 0.9), but particle identification is not 
performed on the other two tracks, at first they are treated as pions. We also require that 
one 7 r° is reconstructed with —5 < 5 77 < 5 from the photons in the signal hemisphere. 
Without any PID for pion candidate tracks on the signal side, possible photon conversions 
can result in a fake event. In order to remove such a contribution, we examine the effective 
mass-squared M 2 e of any two tracks (to which the electron mass is assigned) in the signal 
hemisphere. A sharp peak at M 2 e < 0.01 (GeV/c 2 ) 2 is seen in the Bhabha MC but not in 
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FIG. 2: S^°, S 77 and m 7T + 7T - ri distributions for (a) r~ —> (b) r~ —> /i~rj and (c) t~ —> [i~r( 

decay modes, respectively, prior to final cuts on the signal side (see the text). The data (points 
with error bars), total background MC (open histogram), non-r + r~ background MC (hatched 
histogram), and signal MC (shaded histogram, corresponding to the branching fraction of 5xl0 -7 ) 
are shown. 

the signal and generic t + t~ MC’s, and data also exhibit a tiny peak at M 2 e ~ 0 (GeV/c 2 ) 2 . 
To avoid a large reduction of the signal detection efficiency, we impose a requirement on PID 
in addition to that on M 2 e : events are rejected which have M 2 e < 0.01 GeV/c 2 and V e > 0.6 
for one of the two tracks. As a result, 60% of the Bhabha originated events are rejected 
while the signal efficiency is reduced by 2.3% only. In addition, the following constraints are 
also required: p m i SS > — 1.5m^ iss — 1.0 and p m i SS > 0 . 8 m 2 liss — 0 . 2 , where p m i ss is in GeV/c 
and m miss is in GeV/c 2 . 

The 7r° and rj from 77 are extracted by imposing the conditions —3 < S^ < 3 and 
—5 < S < 3, respectively. The requirements P e / M > 0.9 and p > 0.7 GeV/c are imposed 
on the tracks on the signal side as in the 1-1 prong case. After applying the above conditions 
as well as the requirements on the invariant mass 1.3 GeV/c 2 < M(- m o < 2.0 GeV/c 2 , energy 
difference —1.0 GeV < A E < 0.5 GeV and a loose requirement on 77 ( 7 ') mass m n + 7r - w o < 0.65 
GeV/c 2 < 1.2 GeV/c 2 ), we obtain the resulting number of data and MC events 

summarized in Table an The m^+Tr-rj distribution of the mode with a final state muon after 
these requirements is shown in Fig. [21(c). 

Again, agreement between sideband data and background MC is observed within the 
statistical uncertainty, except for the decay mode r~ —> e~rp 77 —> 7 r + 7 r _ 7 r°. Since in this 
mode PID has better rejection power of hadronic background than in the modes with a final 
state muon, the 7 candidates should be mostly fake ones. In fact, four of the remaining five 
data events are in the 77 mass sideband regions and thus will be rejected at the last stage 
by the tighter 77 and rj' mass requirement 0.5260 GeV/c 2 < < 0.5656 GeV/c 2 and 

0.9027 GeV/c 2 < 777 .^+^-^ < 0.9798 GeV/c 2 , corresponding to the ±3a regions, where the 
resolutions are estimated from signal MC. 


BG ESTIMATIONS AND BRANCHING FRACTIONS 


The signal detection efficiency for each mode is evaluated by using signal MC. We take 
an elliptically shaped signal region in the plane with a signal acceptance, fl, 

of 90%, which is more sensitive than a box shaped signal region. To obtain the highest 





































































sensitivity, all geometrical parameters of the ellipse, such as the position of the center, the 
orientation, and the length of the major and minor axes, are determined to minimize its 
area with the constraint on the acceptance to be 90%. 

As explained in the previous section, an analysis of the background MC distributions 
shows that they agree well with data in the sideband regions. We also confirmed that no 
peaking background structure, which mimics the signal, is observed in the signal region. 

The expected number of background events ( 6 0 ) in the signal elliptical region is estimated 
from sideband data events as follows: for the decay modes t~ —* ^ i~r/ (7 — > 77 ), r~ —* /i”7r° 
and r~ —» , we assume that the background distribution is described by a polynomial 

function in Mi-mo and a Gaussian function in A E over the ± 10 cr region. We determine 
the shape of the function by fitting to the MC events including the blinded region with its 
normalization from the data sideband. Moreover, taking into account the uncertainties from 
assuming the background distribution, we evaluate b 0 for these modes by another method: 
we assume a flat distribution over the ±10cr region of 0 and the regions —0.4 GeV 
< A E < 0.20 GeV, -0.5 GeV < A E < 0.33 GeV and -0.4 GeV < A E < 0.20 GeV for 
the fi ~77 (7 — > 77), and [Erf modes, respectively. The differences between the values 

evaluated by the first and the second method, which are 21%, 10% and 49%, respectively, 
are treated as an additional uncertainty. For the other five decay modes, we evaluate bo from 
sideband data events by simply assuming a flat distribution over ±10cr region of M^- m o-AE 
plane, since the remaining number of events is too small to determine the background shape 
by fitting to those events. The uncertainty on b 0 is quite large (< 100%), particularly in 
the electron modes, because the number of events remaining in the sideband region is very 
small. For t~ —> e~7T°, there are no events in the sidebands, and we set b 0 = O.Ol^g where 
the error is calculated by assuming 2.44 events in the sideband. The resultant b 0 are listed 
in Table EJ 

After opening the blinded region, we find one event in both the r~ —> e~ rj and r~ —> n~r) 
modes, and five events in the r~ —> /i _ 7r° mode, however, no events are found for the other 
modes in the elliptical signal region. Figure U shows the M^-mo-AE plot for the individual 
modes. The resultant numbers of events s are compared to the expected background bo in 
Table HU good agreement can be observed. 


TABLE III: Results of the final event selection for the individual modes: e is the detection efficiency, 
Bmo is the branching fraction for the M° decay, AVA ta and JVM5 are the numbers of events in 
sideband regions for data and MC, respectively, bo is the expected number of background events, 
s is the observed number of signal events and sq is the upper limit on the number of signal events. 


Mode 


Subdecay 

e 

b m ° 

atDATA 

iV side 

* 35 £ 

bo 

S 

so 



mode 

(%) 

(%) 

(ev.) 

(ev.) 

(ev.) 

(ev.) 

(ev.) 

r~ —> e~ 

V 

V 77 

5.68 

39.43 

2 

0 

0.23 ±0.16 

0 

2.3 



r) — > 7T + 7r _ 7r 0 

6.84 

22.6 

1 

0 

0.23 ±0.23 

0 

2.2 

T~ -> 

~V 

V 77 

8.03 

39.43 

9 

9.2 ±2.3 

3.9 ±1.5 

1 

1.4 



77 — > 7T + 7r _ 7r° 

7.15 

22.6 

2 

0.8 ±0.3 

0.60 ± 0.42 

0 

1.9 

r~ —> e~ 

7 r° 

7T° —7 77 

4.70 

98.798 

0 

0.7 ±0.7 

^ 9 
d d 
+ 1 
0 
d> 

0 

2.4 

T~ -> E 

-7T° 

7T° —7 77 

6.36 

98.798 

16 

12.5 ±2.7 

3.0 ±0.9 

5 

6.9 

r~ —> e~ 

rf 

7 / —7 7r + 7T _ 77 

8.51 

17.5 

2 

0.8 ±0.3 

0.28 ± 0.20 

1 

4.2 

T~ -► E 

-rf 

7 / —> 7r + 7r~r/ 

8.41 

17.5 

5 

5.5 ± 1.9 

0.94 ±0.60 

0 

1.6 
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M„, (GeV/c 2 ) 







FIG. 3: The distribution of remaining events over the ±10<r region in the M^- m o — A E plane 
after opening the blinded region. The ellipse is the signal region that has a signal acceptance of 
fl = 90%. The blinded regions ±5<r in M e ~ M o and ±0.5 GeV in A E are indicated by the dotted 
lines. The cross indicates the ±2er M^~ m o and A E ranges. Various symbols show events from 428 
fb -1 of generic t + t~ MC (□), 232 fb" 1 of qq MC with q = u,d,s (A), 160 fb -1 with q = c (o), 
205 fb -1 of two-photon MC (▼), and data (•). 
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TABLE IV: Systematic uncertainties in %. 


Mode r —► 

e rj, 

V -► 77 

e rj, 

rj 

1 ^ 

i 

M V, 
rj —> 37 T 

e 7T° 

1 

O 

e r/' 

ji rj' 

Track recon. 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

M° recon. 

2.0 

4.0 

2.0 

4.0 

2.0 

2.0 

4.0 

4.0 

7r° veto 

5.5 

- 

5.5 

- 

- 

- 

5.5 

5.5 

e ID 

1.0 

1.0 

- 

- 

1.0 

- 

1.0 

- 

H ID 

- 

- 

2.0 

2.0 

- 

2.0 

- 

2.0 

Trigger 

0.5 

0.1 

0.2 

0.1 

0.7 

0.4 

0.1 

0.1 

Beam background 

2.3 

2.1 

2.3 

2.1 

2.3 

2.3 

2.1 

2.1 

Luminosity 

1.4 

1.4 

1.4 

1.4 

1.4 

1.4 

1.4 

1.4 

B m ° 

0.7 

1.8 

0.7 

1.8 

- 

- 

3.4 

3.4 

MC stat. 

1.4 

1.7 

1.1 

1.6 

0.9 

0.8 

1.2 

1.1 

MC models 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

Total 

7.0 

5.8 

7.2 

6.0 

4.2 

4.5 

8.4 

8.6 


Following the Feldman-Cousins method H3, we calculate the upper limit s 0 on the num¬ 
ber of signal events at the 90% confidence level (C.L.) for all modes, as listed in Table Hill 
These values take only statistical errors into account. 

The branching fraction at the 90% C.L. is obtained from the following formula, 


b(t~ -► r M°) = 


s o 


2 e B m °N t 


( 1 ) 


where e is the detection efficiency for the individual modes, Bm° is the branching fraction 
for the M° decay, and N TT = 137.2 x 10 6 is the number of produced r-pairs at 153.8 
fb -1 luminosity, with the cross section of e + e“ —> t + t~ at the Y(45) resonance a TT = 
0.892±0.002 nb calculated by KKMC fisl ]. The values of e and Bm° for each mode are listed 
in Table ITTT1 Systematic uncertainties due to the background estimate bo, and uncertainties 
on 2 eBM 0 N TT , are taken into account by inflating the value of s 0 , as discussed below. 

The systematic uncertainties on the detection sensitivity, 2eBM°N TT , arise from the track 
reconstruction efficiency ( 1 . 0 % for each track) of the tag side track and the signal side 
lepton; 7T°, rj and rf reconstruction efficiency which includes the uncertainties in the tracking 
efficiency for charged pions ( 2 . 0 %) and the n° or rj reconstruction efficiency ( 2 . 0 %); n° veto 
for t~ —> i~r/ and £~rj' modes (5.5%); electron or muon identification efficiency (1.0% for 
electron, 2 . 0 % for muon); trigger efficiency ( 0 . 1 - 0 .7%); beam background effect ( 2 . 1 % for 
1-1 pron g, 2 .3% for 1-3 prong events); luminosity (1.4%); branching fraction of 7T°, rj, and 
rj' (Ref. jlH); signal MC statistics ( 0 . 8 - 1 .7%) and signal MC models (0.5%). Table |IV] 
lists these separate contributions as well as the resulting systematic uncertainties obtained 
by adding all the components in quadrature. The dominant contributions, 7T°, rj and rj' 
reconstruction efficiencies, 7r° veto efficiency and beam background effect, are estimated as 
follows. The contribution of reconstruction efficiencies for pseudoscalar mesons is evaluated 
from the efficiency ratios for data and MC samples using 77 —*• 77 and rj —> 37 t° decays. 
The 7r° veto contribution is also evaluated by the efficiency ratio of data and MC, where 
the 7r° veto efficiency is estimated from the difference of rj reconstruction efficiencies with 
and without the 7r° veto. The beam background effect is estimated from data that is taken 
by the random trigger at the same time as the normal data taking. We initially assumed a 
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uniform angular distribution of r decays for the signal MC. Its possible nonuniformity was 
taken into account by comparing the uniform case with those assuming V — A and V + A 
interactions |25|, which result in maximum possible variations. This effect contributes 0.5% 
shown in the “MC models” line of Table IIVI 

The treatment of systematic uncertainties depends on the estimated background b 0 . In 
cases where b 0 < 1 . 0 , except for r“ —> /i~V, we set b 0 to 0.0 ± 0.0 to obtain a conservative 
upper limit. We use the POLE program |26|) to recalculate s 0 including systematic uncer¬ 
tainties on both &o and the detection sensitivity, assuming Gaussian errors [271. The upper 
limits then obtained from equation (HD are summarized in Table El 


TABLE V: Upper limits on branching fractions. 


Mode 

Subdecay mode 

U.L. of 8 @ 90% C.L. 

t~ —> e~rj 

V 77 

4.0 x 10“ 7 


ij —> 7r + 7r _ 7r° 

5.8 x 10 -7 

t~ —> e~r] 

combined 

2.4 x 10 -7 

t~ -»• n~r] 

V * 77 

2.3 x 10" 7 


77 — > 7r + 7T~7T 0 

5.5 x 10" 7 

t~ -»• PW 

combined 

1.5 x 10“ 7 

T~ e _ 7T° 

7T 0 —> 77 

1.9 x 10“ 7 

T~ -4 p“7T° 

7T 0 —> 77 

4.1 x 10“ 7 

t~ —> e~r( 

rj —>• TT + TT~r] 

10 x 10- 7 

T~ -»• n~rf 

rj —>• TT +/ K~r] 

4.7 x 10“ 7 


DISCUSSION 


One can see from Table El that for the r —» e r/' and r —> /i rj' modes, the resulting 
90% C.L. upper limits are B(t~ —> e~rj') < 10 x 10 ~ 7 and B{t~ —> /i - ?/) < 4.7 x 10 -7 ; 
these are the first experimental limits on these modes. For the r~ —» e~r] and —> p - r/ 
modes, the results from two different final states, 77 —> 77 and 77 —> 7 r + 7 r _ 7 r°, are combined. 
The resulting 90% C.L. upper limits for these two modes and t~ —> e~7T° and r~ —* /i _ 7 r° 
are B{t~ —> e~r /) < 2.4 x 10 -7 , B(r~ —> /U 77 ) < 1.5 x 10” 7 , B(t~ —> e _ 7 r°) < 1.9 x 10 - ' 
and B(r~ —> /i _ 7 r°) < 4.1 x 10“' and improve upon the previous CLEO measurements by a 
factor of 34, 64, 20 and 10, respectively. 

We can restrict the allowed parameter space for tan/3, using the relationship derived 
by M. Slier jd] in a seesaw MSSM with a specific mass texture: 


B(t~^> fTv) = 0.84 x 10“ b x 


/tan/A 6 /100 GeVV 
) V m A ) 


( 2 ) 


where mA is the pseudoscalar Higgs mass and tan (3 is the ratio of the vacuum expectation 
values ((H u )/(Hd)), as shown in Fig. Q] where our boundary is indicated for the 90% C.L 
by the shaded region. Figure 01 also shows the 95% C.L. constraints from our experiment 
with a dashed line and high energy collider experiments at LEP [20] and CDF |29|. Our 


12 














experiment has a sensitivity competitive with that of the CDF experiment which searched 
for pp —> A/(pbb —> bbbb reactions, where 0 is a CP-even neutral Higgs state and A is a 
CP-odd state in the MSSM. 

The improved sensitivity to rare r lepton decays achieved in this work can also be used 
to constrain the parameters of models with heavy Dirac neutrinos j@, 6]. In these models the 
expected branching ratios of various LFV decays are evaluated in terms of combinations of 
the model parameters. These combinations, denoted y Te and y T/1 for r decays involving an 
electron and a muon, respectively, can vary from 0 to 1. Our r~ —> e“7r° result sets a 90% 
C.L. upper limit y Te < 0.26, the most restrictive bound on this quantity. The corresponding 
limit from our r~ —> p _ 7r° result is y T/1 < 0.87: a somewhat better bound y Tfl < 0.65 can 
be set from the r~ —> /i + decay based on the experimental limits from BaBar dOf and 

Belle 0. 

In summary, using a data sample of 153.8 fb _1 collected with the Belle detector we ob¬ 
tained new upper limits on the branching fractions of semileptonic LFV r~ decays involving 
pseudoscalar mesons 7r°, y and y'. They range from 1.5 x 10~ 7 to 10 x 10“' for the six decay 
modes studied and are 10 to 64 times more restrictive than previous limits. 



m A (GeV/c 2 ) 


FIG. 4: Experimentally excluded mA — tan/3 parameter space. The 90% C.L. result of this exper¬ 
iment on £>(t~ —> y~y) using the relation @ is indicated by the shaded region together with the 
95% C.L. regions excluded by this experiment (dashed line), LEP 0 and the CDF 24, 0- 
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